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The partition function of two dimensional massless staggered fermions interacting with U(N) gauge fields is 
rewritten in terms of loop variables in the strong coupling limit. We use this representation of the theory to devise 
a non-local Metropolis algorithm to calculate the chiral susceptibility. For small lattices our algorithm reproduces 
exact results quite accurately. Applying this algorithm to large volumes yields rather surprising results. In 
particular we find m, 7^ for all N and it increases with N. Since the talk was presented we have found reasons 
to believe that our algorithm breaks down for large volumes questioning the validity of our results. 



1. INTRODUCTION 

Understanding the chiral limit of lattice gauge 
theories with fermions still presents a formidable 
challenge for numerical work. Conventional al- 
gorithms like the hybrid Monte-Carlo or multi- 
boson techniques encounter difficulties in this 
limit. In this article we discuss how a new loop 
representation helps one to solve the chiral limit 
of strongly coupled U(N) gauge theories interact- 
ing with staggered fermions in two dimensions. 

The chiral limit of strongly coupled gauge the- 
ories involving staggered fermions in two dimen- 
sions contain some unresolved questions. It is well 
known that staggered fermions preserve a U(l) 
chiral symmetry at finite lattice spacings. Using 
large N (color) Q and large d (dimension) || 
approximations, it is possible to show that this 
chiral symmetry breaks spontaneously at strong 
couplings. There are also some rigorous results in 
four dimensions jjj. As a consequence, at least in 
d > 3 the model is expected to contain a single 
Goldstone boson (pion) which is exactly massless 
in the chiral limit. What happens in two dimen- 
sions? Numerical simulations suggest the exis- 
tence of a massless pion [0. This is somewhat 
surprising since in d — 2 the Mermin- Wagner the- 
orem forbids the breaking of a continuous sym- 
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metry and the pion cannot be a Goldstone bo- 
son. On the other hand in two dimensions we 
know that a U(l) symmetry is special. Kosterlitz 
and Thouless have shown that a theory with this 
symmetry can contain both a massive phase and 
a massless phase separated by a transition. So an 
interesting question is whether the strongly cou- 
pled U(N) lattice gauge theory with staggered 
fermions contains a massless pion in the chiral 
limit or not and if it does is the long distance 
physics in the same universality of the 2d X-Y 
model. 

2. LOOP REPRESENTATION 

About two decades ago, it was shown that a 
U(N) gauge theory involving a single flavor of 
staggered fermions in the strong coupling limit in 
any dimension is equivalent to a statistical me- 
chanics of monomers and dimcrs @. The parti- 
tion function is given by 

[n,b] x,f* x 

where x represents a lattice point in the d di- 
mensional lattice, n = 1,2, ...d the direction, p, 
the unit vector in the corresponding direction, 
n x = 0, 1, ...,N is the number of monomers lo- 
cated on the site x and b Xjf ^ = 0, 1, 2, N is the 
number of dimers located on the bond connecting 
the site x and the nearest neighbor site x + jl in 
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the /i direction. We will assume b x>ll = b x +p, t -n 
for simplicity. The configurations [n, b] that con- 
tribute to the partition function are constrained 
by the relation 

n x +Y^ &*,|. + &x,-m = N (2) 

on each site x. 

In the chiral limit n x is forced to be zero and it 
is difficult to design a local algorithm for the par- 
tition function that obeys the constraint. Perhaps 
for this reason the chiral limit remained unex- 
plored. In order to avoid this problem we rewrite 
the partition function given in ([!]) in terms of new 
"loop" variables. For the case of N = 1 this was 
recently explained in Q. We first shade plaque- 
ttes of the two dimensional lattice as in a chess 
board. We then introduce new (dashed) dimers 
d x ,fi = 0, 1, 2..JV on each bond such that each 
corner x of any shaded plaquette satisfies 

b x>fi + d Xtfl + b x%v + d x . v = N (3) 

where x + ft and x + v are the two nearest neigh- 
bor sites of x also belonging to the same shaded 
plaquette. In (||) ft and v can be negative if the 
direction of the nearest neighbor site is backward. 
Of course there are many ways to introduce the 
d X)A( bonds that satisfy (||). It is easy to check that 
for a shaded plaquette with bonds 61,62,63,64 
in a cyclic order, the total number of ways one 
can introduce the dashed bonds obeying the con- 
straint @ is (7V+l-max(61,63)-max(62,64)). 
Thus in order to keep the partition function the 
same we re-weight each extended configuration 
with l/(N + 1 - max(61, 63) - max(62, 64)). In 
this construction it is easy to see that each site 
touches N — n x solid (original) dimers and N + n x 
dashed dimers. When n x — one can at ran- 
dom connect every solid dimer to a unique dashed 
dimers in Nl ways. In the chiral limit this re- 
sults in loops made up of an alternate sequence 
of solid and dashed dimers. In figure [l] an N = 2 
configuration of solid dimers on a 2 x 2 lattice is 
extended to a loop configuration using the above 
rules. An update of the system involves flipping 
of a loop which means converting all solid dimers 
into dashed dimers within the loop and vice versa. 
Typically, a point is selected at random and the 
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Figure 1. An N — 2 loop configuration on a 2 x 
2 periodic lattice. The two plaquettes are shown 
separately for convenience. The numbers at the 
end of the dimers indicate connections. There are 
two loops in this configuration. 



loop connected to it is flipped using a Metropolis 
decision. Since there are (b + d)\ / (d\b\) number of 
equivalent dimer configurations with 6 solid and 
d dashed dimers on a given bond, it is necessary 
to take this factor into account in the Metropolis 
decision to ensure detailed balance. 

3. CHIRAL SUSCEPTIBILITY 

To address the question of the existence of a 
massless pion in the chiral limit one can measure 
the chiral susceptibility x as a function of the 
lattice size L. Since it is the first derivative of 
the chiral condensate with respect to the fermion 
mass it should diverge as a power of L if there 
is a massless pion in the theory. For example if 
the chiral symmetry is spontaneously broken we 
must have x oc L d where d is the dimensions of 
the lattice. In the case of two dimensions where 
we expect chiral symmetry to remain unbroken, 
a massless pion leads to a divergence of the form 
X oc Z/2 — r/) where the KT prediction is that 
< Tj < 0.25. Another example is the continuum 
two flavor Schwinger model where rj = 1. 

To measure the chiral susceptibility, we intro- 



duce two monomers into each loop configuration 
in all possible ways. Since the loops are made 
up of an alternate sequence of solid and dashed 
bonds, introducing a monomer introduces a de- 
fect in the pattern. However, a single defect is im- 
possible in a loop. Thus, in a two monomer con- 
figuration both monomers must be on the same 
loop. It is easy to see that a loop with two 
monomers can be generated by partially flipping 
a loop. Using this observation we find the chiral 
susceptibility by summing the weights of all the 
new configurations generated during the loop up- 
date. We have checked that our measurement of 
the susceptibility produces accurate answers on 
small lattices. Table 1, compares our results to 
exact calculations. 

Table 1 

Susceptibility: algorithm vs. exact results. 



N 


Lattice Size 


Exact 


Algorithm 


1 


8x8 


5.27221... 


5.2722(2) 


2 


4x4 


6.40205... 


6.402(2) 


3 


4x4 


14.1595... 


14.159(8) 


5 


2x2 


9.57386... 


9.574(7) 


30 


2x2 


338.534... 


338.2(8) 



When we apply this algorithm to larger lattices 
we find that small clusters often get updated and 
occasionally even loops of large sizes get flipped. 
Our results for the chiral susceptibility for var- 
ious lattice sizes is plotted in figure ^. We can 
draw three conclusions from the plot. First, as a 
function of L there is no evidence of a divergence 
for any AT. Thus the pion is always massive. Sec- 
ond, since the value of L where the susceptibility 
begins to level gives a rough indication of the in- 
verse of the pion mass, the pion mass increases 
with AT. Finally, the value of the chiral suscepti- 
bility diverges with N, signaling the formation of 
a condensate. However, our results suggest that 
this is not driven by long distance fluctuations. 

Since the talk was given we have continued to 
study our algorithm more carefully. It now ap- 
pears to us that there are configurations which 
we sample exponentially rarely as the volume in- 
creases but whose contribution is still finite to 
the final result. This problem begins to show up 
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Figure 2. Chiral Susceptibility as a function of 
the lattice size L for various values of N. 



mildly around L — 32 for N — 1 through very 
large fluctuations. Of course we had collected 
very very large statistics to be able to control our 
error bars. However, we think the error analysis is 
not reliable since we may have uncontrolled auto- 
correlation times. For this reason we believe the 
error bars in figure || and in the graphs presented 
in H are incorrect. This also invalidates the con- 
clusions made above and in |J. Thus, although 
our representation is a useful first step and per- 
haps gives a qualitative picture for small volumes 
more work is needed before we can answer ques- 
tions we raised in the first section. 
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